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The magnetic properties of a charged perfect quantum gas in the independent electron approxi- 
i-C . mation and confined to a box A have been extensively studied in the literature. One of the central 

problems has been to establish the thermodynamic limit for the magnetization and magnetic 
. susceptibility, see e.g. [T], [2], [3J, [3] ,0, [BJ, [5] and references therein. 

Briefly, the technical question is whether the thermodynamic limit (A — > M 3 ) commutes with 
the derivatives of the grand-canonical pressure with respect to the external constant magnetic held 
B. 

A general way of proving this thermodynamic limit has been already announced in [5] and 
[BJ, and the main ingredient consisted in applying the magnetic perturbation theory to a certain 
Gibbs semigroup. The strategy of the proof, which works not only for the first two derivatives, 
but also for derivatives of all orders, were outlined in [3J. In [3] we proved in detail the pointwise 
thermodynamic limit near z = 0. This paper is the last one of the series, it contains the complete 
proof of the uniform bounds on compacts needed in order to apply the Vitali Convergence Theorem 
(N : (see [TUJ). 

Now let us formulate the mathematical problem. The box which contains the quantum gas 
will be the cube A C M 3 of side length L > centered at 0. The constant magnetic field is 
B = (0, 0, B), with B > 0, oriented parallel to the third component of the canonical basis in K 3 . 

We associate to B the magnetic vector potential Ba(x) = x%, xi, 0) and the cyclotronic 
frequency uo = -B. In the rest of the paper, u> will be a real parameter. The one particle 
Hamiltonian we consider is the self-adjoint operator densely defined in L 2 (A): 

H L {u):= l -{-iV -^a) 2 , (1.1) 

corresponding to Dirichlet boundary conditions. 

One denotes by -Bi(L 2 (A)) the Banach space of trace class operators. At oj > fixed, the 
magnetic Schrodinger operator Hl{u) generates a Gibbs semigroup {Wl(0, w)}p>o where: 

W L (P,u;):=e-^ u \ \\W L (f3,Lo)\\ Bl < — ^- T , (3 > 0. (1.2) 
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Let us now introduce the grand-canonical formalism. Let = l/(kT) > be the inverse 
temperature, fj, € K the chemical potential and z — the fugacity. Let K be a compact 
included in the domain D + := C\ [e~,oo[ for the Bose statistics and £>_ := C\] — oo,—e~] for 
the Fermi case. 

Fix wo > and a compact real interval £1 containing ojq. For a fixed f3 > 0, one can find a 
simple, positively oriented, closed contour Ck C D± whose interior does not contain 1 in the Bose 
case or —1 in the Fermi case, and such that 

sup sup sup sup || [£ — zW/l^w)] -1 !] = M < oo. (1.3) 

L>1 wGO ££C k zeif 

Details may be found in [5] for the Bose case, but the main idea is that the spectrum of Wl((3,u)) 
is always contained in the interval [0, e - ' 9 "/ 2 ] and one can apply the spectral theorem. 

We then can express the grand canonical pressure at wo as follows (see e.g. [5] for the Bose 
case): 

P L (f3, z ,u») = — € — f dg ln(1 ~ ^ Tr [(£ - zW L (0 ) w o ))- 1 zW L (/3,u o j\ ■ (1.4) 
ziirpij j Ck ^ 

where e = 1 for the Bose gas, and e = — 1 for the Fermi gas. 
For uj E R, £ £ Ck and z £ K, introduce the operator: 

g L {p,z,Z,u) := [Z-zW L {p,u)\- l zW L {fl,u). (1.5) 

This is a trace class operator which obeys (use (|1.3[) and (|1.2[l ): 

\\gL<fi,u,M\\ Bl < (su P |z|)— — -. (1.6) 
zeK (27rp)2 

uniformly in w £ fi. Because w — > Wl(P,u>) is a Bi-entire operator valued function in u> (this 
result was first obtained in [TJ and then refined in [3]), then using (|1 .3[) one easily shows that the 
map 

]0, oo[9 uj I— > Tr ,gi(/3, -z, £, cj) € C 

is smooth, with derivatives which are uniformly bounded in £ and z. Thus for every N > 1 and 
zgifwe can define the generalized susceptibilities at cjq by 



X&(A^o) := ^ ~„ (P,z,uj ) (1.7) 



dt, ~, a , n (A^C^o). 



2inf3L 3 J Ck s £ 0a/ 

From this discussion one can see that the pressure as well \l (A "j w o)j > 1 are analytic functions 
on D + (or D_). 

Now let us describe the case when L = oo. The thermodynamic limit of the pressure exists 
and is uniform on compacts like K. If uj > 0, its actual value is (see [2] ): 

Poofo *,<*>) :=o;q — * £/| /2 , (L8) 

where /^(C) are the usual Bose (or Fermi) functions for e = 1 (or e = —1): 

£(C) : = -fK / * -i 7 _ t , (1-9) 

analytic in C \ [1, oo[ (or C\] — oo, —1] ) if e = 1 (or e = —1). If \Q\ < 1, they are given by the 
following expansion: 

xi ,n— l/-n 

£(0 = £ ! -^-. 

n=l 
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Now it is rather easy to verify that for any N > 0, the multiple derivative Pqo(/3, •, Wo) exists 
and defines an analytic function on D + (or D ). The main result of [3] established the pointwise 
convergence, 

lim dZP L ((3,;u >0 ) = d»P 00 ((3,;u >0 ):=xZ(P,;Uo); \z\<l. 

L — >oo 

Remember that we want to apply the Vitali Convergence Theorem (see [10] or [3]). Therefore, 
in order to conclude that Xl (A z i w o) converges uniformly to Xoo(P> z i w o) for au z ^ K, the only 
remaining point is to get the uniform boundedness w.r.t. L. More precisely, we will prove: 

Theorem 1.1. For all N > 1, for all (5 > and for all uj > 0, 

sup sup \xl (P,z,u))\ < comt(0,K,w,N). (1.10) 

L>1 zeK 

Then putting this together with the pointwise convergence result near z = of [4], the final 
conclusion would be: 

lim sup (/?,«. Wb) I =0. (1.11) 



L- 



: z£K 



Remark 1.2. Having uniform convergence (|1.11|) «ra£/i respect to z allows us to prove existence 
of the thermodynamic limit for canonical susceptibilities (see J3[ \bT/) 

This paper is devoted to the proof of Theorem 11.11 Note that the theorem is an immediate 
consequence of the following estimate: 



sup sup 

£eC K z£K 



d»Trg L 



< L 3 const((3,K,N,L0 ), (1.12) 



doj N 

which would imply via (11. 7|) that the generalized susceptibilities are uniformly bounded in L. 
1.1 Strategy of the proof 

From now on, we omit the parameters £ and z in the definition of g^ in order to simplify notation. 
Fix (3 > and ujq > 0. Let Q C K be a compact interval containing loq. If w € f2, we denote by 
<5u := a; — o^o- The main idea of the proof is to derive an equality of the following type: 

N 

Tr 5L (/3,^) = Tr 5i (/?,wo)+^(M j ^(Awo) + (Su) N+1 K L (P, w, N), (1.13) 

i=i 

where the coefficients aj(/3,uJo) grow at most like L 3 uniformly in £ and z, while the remainder 
1Zl(P, N) is a smooth function near u>q. Then since we know that Tr Ql(P, ■) is smooth, we must 
have 

d N Tr g L 
dujN (A Wo) = N\a N {f3,uj ), 

and this would finish the proof. In order to achieve this program, we will have to do two things. 

First step: with the help of magnetic perturbation theory we will find a regularized expansion 
in Su) for of the form 

N 

g L (fi, u) = J2( S ") n 9L,n(P, w) + Rl,n((3, u, N), (1.14) 

which holds in the sense of trace class operators, and the remainder has the property that 
(5u)x+iRl,n((3, uj) is smooth near lo in the trace class topology. The operator-coefficients gL, n 
will still depend on to, but in a more convenient way. That is, they are sums, products, or integrals 
of products of regularized operators, see ()2.2|) . This result is precisely stated in Theorem [ 
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Second step: show that for each < n < N we can write 

N 

Tr 9lAM = YJ&tf 8L *>»(P><*i) + (Su) n+1 K l , n (I3,lj,N), 

3=0 



(1.15) 



where the remainder 1ZL, n {(3, ■) is smooth near wq. Now the coefficients S£j, n (/3, ljq) are finally 
independent of u>, and grow at most like L 3 . This is done in the last section. 
Finally, if we combine (|1.15l) with (|1.14p , we immediately obtain (|1.13p . 

Now let us discuss why a more direct approach only based on trace norm estimates cannot 
work. Recall that the map to — > Wl(/3,u>) € B\ is real analytic, hence d a W N L is well defined in 
Bi(L 2 (A)), and we have the estimate (see [T] and [I]) : 



1 d Wl 
iV! 8uj n 



L 3+N {l + P) sN 



Bi 



[3 



2 rJV-1" 



(1.16) 



where cn is a positive constant which depends on N, u>q and s. Now if we use the Leibniz rule of 
differentiation for the product which defines the operator gL(f3,ui) (see (|1.5j) ). and estimate traces 
by trace norms we obtain: 



d N Tv g L 



dLO N 



08, wo) 



< 



d N 9L 



8uj n 



(13, wo) 



< L 



3+N 



const (/3, if, j N). 



(1.17) 



Now this is definitely not good enough, and we have to find a more convenient expansion, as 
described in (|1.15[) . This will be done in the next sections. 



2 Regularized expansion of Wl 

It has been shown in [3] and [3] that by using gauge invariance one can control the linear growth 
of the magnetic vector potential a. The price one pays is the introduction of an antisymmetric 
phase factor, which disappears though when one takes the trace. Let us now show how this works 
for the operator Wl- 

Fix luq > and f3 > 0. Let u> € C and 5u> as above. Let us define the magnetic phase: 

<Kx,x') :=x-a(x') = -(x 2 x[ - Xl x' 2 ) = -</>(x',x); (x,x')eA 2 . (2.1) 

If r(cjo) is a bounded operator with an integral kernel £(•, •, coq), then the notation T(u>) will refer 
to the regularized operator associated to T(u>o) which has the kernel: 

i(x, x', u) := e^( x ' x ')t(x, x', u ); (x, x') e A 2 . (2.2) 

We will very often use the Schur-Holmgren criterion of boundedness for integral operators (see 
[H]), which states that if T has an integral kernel t(x, x'), then: 




We denote by Gl(-, •, (3, lu) the kernel of Wl((3, uj). We define two other bounded operators 
R\ t L and i?2,L by their kernels, 

i?i, L (x, x', (3) := a(x - x') • [iV x + oj a a (x)] G L (x, x', (3, ui ), 

fl 2 , L (x,x',/3) := a2 ( x ~ X,) C? L (x,x / , j 9,a;o); (x,x')eA 2 . (2.4) 
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Then consider the corresponding regularized operators Wl, R\,l, Ri.l- Let us state here two 
important estimates, the first one is just the diamagnetic inequality, while the second one was 
obtained in [5], 

|G L (x,x',/3, Wo )| < G oo (x,x , ,/?,0) = (2^)3/2 ex P (~ |X 2/? X/| ) ' (2 ' 5) 

|[zV x + ^ a(x)]G L (x,x',/3,^ )| <C(1+lu ) 3 { G OQ (x, x', 8/3, 0), (2.6) 

on A 2 , where C > is a numerical constant. A straightforward application of the Schur-Holmgren 
criterion gives us the following operator norm estimates 

||WiiG9,w)||<l, \\Ri tL (0,uj)\\ < Co(l + ^ ) 3 (l + /3) 5 , (2.7) 

where i = 1,2 and Co > is a numerical constant. 
For ix, i n £ {1, 2}, define 

Dn{P) = {(ri, ...,r„) 6l° : < r„ < ... < n < /?}. (2.8) 

Introduce the operator norm convergent Bochner integrals: 

I n ,L(h, ...,£„) (/?,u;) := / drWi^-rijwjEij.i^i -T 2 ,a;) 

• Ri2,l(t 2 - r 3! o;)...Ei n _ 1 ,i(T n _i - r„,cj)^i n) L(r„,w), (2.9) 

and 

J ntL {(3,uj) := / drW L {(3 - t 1 ,u)R l (ti - r 2 ,w) 

• -Rl(t 2 - r 3 ,w)...iii(T n _i - r ni uj)R L (T n ,uj). (2-10) 

Here we used the notation: 

Rl(P,v) = (5u)R ltL (p,w) + (M 2 ^2,l(/3,w). (2.11) 

By (|2.7|) the operators I n ,L(ii, •••,*«.) and J nj L belong to B(L 2 (K)). We will show below that in 
fact their belong to I?i(L 2 (A)) and their trace norm is of order L 3 . Denote by Xni^ii ■■■■> in) the 
characteristic function of the set 



{(*!,...,.«)€ {1,2}" : J2% k =j}. 



fe=l 



Proposition 2.1. Fix cl>o > 0, and iV > 1. Set <5w = u> — wo, w S C. Then we have the following 
identity in Bi(L 2 (A)): 



N 



W l (P,uj) = W L (M +^(M"^(i3,w) + 4 1) JV (/3,w), (2.12) 



n=l 



ere — > Tg^prn -^i, jy(A CJ ) * s a smooth B\ operator valued function in u>. 
The coefficients of the above expansion are given by: 

n 

W L , n (M :=^(-l) fe X^{ii,---,ik)hAh,...,ik){P,u): (2.13) 

fc=l ijG{l,2} 

and the reminder reads as: 

2N N 

R%{[3^):= (M n £("!)* £ Xt{ix,...,ik) (2.14) 

n=W+l fc=l i s e{l,2} 

■ hAh, -,ik){P,u) + {-l) N+1 J N+hL ([3,Lo). 
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Proof. This proposition was proved in [2] in the sense of bounded operators. Again by using 
Remark 3.4 in |4|, the operator Wl((3, oj) as well as Wl,t»(A w ) belong to Bi(L 2 (A)). The same 
argument holds for the remainder R^ N {f3,uj). The fact that the operator uj — > Tg^pm R^' N {[3,w) 
is a smooth i?i operator valued function near ujq follows from definitions (|2.14ll . (|2.9|) and ()2.10p . 

□ 

Remark 2.2. If is important to note that the coefficients Wl,n in (MUM) s ^ depend on uj, but 
only through the magnetic phases. Let ui G R. From the estimates (|2.7[) and the definitions in 
(|2.9p and ()2.10p . and after integration over the r variables, clearly we get as bounded operators: 

\\R^ N (f3,oj)\\ < const \5u\ N+1 , (2.15) 
where the constant is uniform in L > 1. 



3 Regularized expansion of gi 

We will now try to get a similar expansion for g^, aiming to obtain (11.14jl . Fix [3 > and wo > 0. 
The parameters z and £ which enter the definition of g^ are fixed as in (11. 5p . and the estimates we 
make must be uniform w.r.t. them. Here uj is real and Suj is as above. It was shown in [5] that: 

Lemma 3.1. Let uj G K. The trace class operator g^ admits a continuous integral kernel. More- 
over, there exist two positive constants C and a, both independent of L, such that the integral 
kernel satisfies 

| 9L (x,x'; W )|<Ce-l x - x 'l. (3.1) 

Looking at (| 1 . 5[> . we see that we need a regularized expansion for the operator (£— zWl((3, 
Using (j2~T5"l) with N = 1 we get for uj G K, 

\\W l {(3,uj)-W l {P,uj)\\ <Cx \8uj\, (3.2) 
where Ci is a L independent constant . Hence choose \8uj\ is small enough such that 

dM\Suj\ (sup \z\) < 1/2, (3.3) 

z£K 

C\ and M being respectively defined in (|3.3p and (II. 3p . Then the operator 1 — z(Wz, — — 
zWl)' 1 is invertible and its inverse has a norm less than 2. Hence by choosing fl to be a small 
enough interval around ujq we get 



sup sup sup sup 



[C-zW^/?,^)]- 1 <2M. (3.4) 



Then we can write: 

ii-zWdP^))- 1 (3.5) 

OO 

= (£ - zW L (P, uj))- 1 ]T z"{[W L (/3, W ) _ w L (P, uj)} (£ - zWU/3, uj))- 1 } 71 

71=0 



and thus we reduced the problem to the study of (£ — zWl((3, uj)) 1 . 

In order to get a convenient expansion for this inverse, we need to introduce some new notation. 
Let x, x' G A. Define for TV > 1: 

r L!A r(x,x ;(J) := -z — G L (x, y; uj )g L (y, x ; uj )dy, (3.6) 
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where 

fl(x, y, x') := 0(x, y) + </>(y, x') + 0(x', x) = ie 3 ■ [(y - x) A (y - x')]; (3.7) 

with e3 = (0, 0, 1), denotes the magnetic flux through the triangle defined by x, y, and x'. Similarly 
let w € K and define the bounded operator rx(/3, c<;) and fx(/3, u)) whose kernel is given by: 

r L (x,x';/?,a,) = -z J ( e *<"fl(W) - l) Gx(x, y; f3, wo)ffx(y, x'; (3, LUo)dy, 

fi(x, x'; 0, u) = e^ x ' x 'V L (x, x'; (3, u) (3.8) 

Notice that fx does not coincide with the regularization of rx given by ()2.2|) . The operator rx 
and fx are related to the operators rx,jv and ?l,n respectively by 

oo oo 

r L (f3,w) = ^(fe) fe rx, fc (/3); fx(M = ^jto) k f L M- (3-9) 

k=l k=l 

Note that by using the Schur Holmgren criterion, the Lemma l3~l"1 (|2.5|) and the fact that fl(x, y, x') 
is bounded from above by L? on A, we have the estimate: 



N oo 

r L -^(MV.fc = (<M fc?, x,fc 

fe=l fc=iV+l 



< const e L \5u\" + \ (3.10) 



for some numerical positive constant. The same estimate again holds true for fx(/3, cj) and the 
corresponding series given in (|3.9|) . 

Let us now give some more precise estimates on the norms of these operators. 

Proposition 3.2. Fix N > 1. There exist a positive constants C<x independent of L > 1 swc/i £/ia£ 
/or a// weR, 

max|max||rx, fe (/3)||, (M^IM/?, < C 2) (3.11) 

and 

max|max||rx, fc (/3)|| B2 , IM^IMA w)|| Ba | = C 2 • L 3 / 2 . (3.12) 

These estimates also hold true for the regularized operators in the sense of (|2.2[) and fx(/3, w). 

Proof. First, note that |fl(x, y,x')| < |x — y| |y — x'|, see (|3.7j) . The kernels present in the y 
integral are localized near their diagonal, see (|2.5p and (|3.1|l . By extending the integral with 
respect to y over the whole K , then using a fraction of the exponential decay in order to bound 
the polynomial growth from the flux, we obtain a constant independent of L such that 

|rx, fc (x,x';/3)| < Ca-e-T^'l; (x,x')eA 2 , 1 < k < N. (3.13) 

The same estimate holds for rx- Now we can apply the Schur-Holmgren criterion (|2.3p and get 
(|3.1ip . The Hilbert-Schmidt estimates is also straightforward. 

□ 

The next proposition gives the necessary expansion of (£ — zWl([3,uj)) . 

Proposition 3.3. Fix N > 1 and loq > 0. Then if \Su\ is small enough, the following identity 
holds in B(L 2 {A)): 

tt-zWLi^u))- 1 =r 1 (l + 5L(/3,^)) [l + r^CS.w)]- 1 (3.14) 
JV 
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where Sl,n is given by 

N 

:= r 1 E X^i,--,y(l+5i(M) 

n=i (ii,...,i„)e(w*)» 

•rc, il C9,w)...rL liB C9,w) (3.15) 

(2) 

where the remainder R L N (f3, to) has the property that the bounded operator valued function oj — > 

(8oj)~ n R^ N {(3, oj) is smooth around luq and moreover, there exists a constant (possibly) depending 
on L such that 

\\R%{(3,u>)\\ < const \6u\ N+1 . (3.16) 
Proof. We start with the following resolvent equation, 

(t-zWLiP^))- 1 =C 1 + (Z-zW L (l3,u;))- 1 zW L (P,u J )C 1 - (3.17) 

Now the next identity is very important, and it is obtained by a straightforward calculation from 
and the definition of gi, (see also Proposition 13 in [5]) 



K - *W L <fi, u)]g L (p, lo) = zW L (f3, oj) + f L (f3, oj). (3.18) 
If one multiplies with an inverse both sides of the above equality we get: 

(C - zW L (fi, oj)Y x zW l {I3, oj) = g L {(3, w) - (£ - zW L (fi, uj))' 1 ^, oj). (3.19) 
We know from Proposition 13.21 that we can find a constant C2 independent of L such that 

\\r L {P,uj)\\<C 2 \8oj\. (3.20) 
Let us use (|3.19l) in (|3.17[) . and isolate the inverse we are interested in 

K - zWU^oj)}- 1 ^ + r 1 riC8,«)] = r^l + ftLCS.w)]. (3.21) 

Now if 1 5a; I is small enough, 1 + £ —1 r£,(a;) is invertible and (|3.14p follows. Moreover, expressing 
the inverse by a finite Neumann-type expansion, 

JV 

[i + rvr 1 = 2>*)~M + 1 1 + r^-H-o-^^f 

fc=0 

and using (|3.9p we can identify the operators Sl,n{P,w) as given in (|3.15| . while the reminder 
reads as 

R^ N (M := (-0- (Ar+1) K " zWUM]- 1 ?"* 1 ^) +SiC8,w)] 

00 AT 

• £ (&;)*X)(-r 1 ) B £ xS(ii..",^)ri,i 1 09,o;)...r £ , <w 09 > a;). (3.22) 

k=N+l n=l (ii,—,in)e(iV*) n 

Let us now identify the term in (<5a;) Ar+1 which appears in the estimate (|3.16j) . For the first term 
of the remainder it comes from ()3 . 20(1 , while for the second one it comes from the fact that the 
series begin with the index N + 1 (see (|3.10[l ). □ 
We are now ready to give a convenient expansion for the operator [£ — zWl(P, u)]^ 1 . First we 
need some new notation. We introduce the following operators, 

S L ,oG9,w) -^[I + SlCM], (3.23) 

N 

T LiN ((3,ui) := 5^2™ X^+i(«o, ii,«i, -, jn,in)S L ,i (f3,u) 

n=l 0<i k <N,l<j k <N 

•W L &1fi ) w)S LM {fl ) w)...W LJn {fl ) w)S L , in {fl ) w) 1 N>1. (3.24) 

Since the operators Wz,j and defined in Propositions 12. ll and [3~3l arc uniformly bounded in 
L, this is also true for Tl n- 
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Corollary 3.4. Fix N > 1 and ujq > 0. If \6u>\ is small enough, then the following identity holds 
mS(L 2 (A)): 

N 

[Z-zWl&u)]- 1 = [£- zWl&u)]- 1 +J2(Su) n T L , n (l3,u;) + R^&w), (3.25) 

n=l 

where the remainder R L N (/3,ui) has the property that the bounded operator valued function u> — > 

(5u)) N ^LN (A u ) * s smo °th near loq, and there exists a constant (possibly) depending on L such 
that: 

\\RW N (p,u)\\ < const \5w\ N+1 . (3.26) 

Proof. The result follows after inserting the estimates from the previous proposition into formula 
([3~5]) . having used the notation introduced in ([333"]) . (|3~2"4"|) . ([2~T2"f and PT5| . The rest is just a 
tedious bookkeeping of various terms. □ 
We finally are in the position of writing "the right" expansion for the operator <?i(/3, u>) as 
announced in (|1.14p . 

Theorem 3.5. Fix N > 1 and loq > 0. // \ 6u>\ is small enough, then the following equality takes 
place in B 1 (L 2 (A)): 

N 

g L (M = g L fi{P^) + Y.^T 9lAP^) + Rl ] n{P^), (3-27) 

n=l 

where 

g L ,oW,o>) := [e-zW^CMr^iCM- (3-28) 
and are given by (N > 1), 

N 

g L , N (P,w) := J2 [S L . N -n(l3,uj)zW L , n {(3,oj) +T L ^u)zW L ,N- n {M\ , (3-29) 

n=l 

where Wl.o '■= Wl and the remainder (g^w Rl^n(P> ^) has the property that the B\ operator valued 
function 

^ — * (5 \ N jv ^ L? ) ^ STiiooth near ujq and there exists a positive constant (possibly) 
depending on L such that: 

\\R^ n (P,lj)\\ Bi < const \Slo\ n+1 L 3 . (3.30) 

Proof. First we multiply the Bi(L 2 (K)) expansion (|2.12|) of the semigroup with the expansion 
(f3~2"5]) of the resolvent valid in B(L 2 (A)). Thus one obtains in Bi(L 2 (A)), 

N N N 

+ ^(<5w)" T L , n (0, uj)zW l ((3, u) + J2 E( fe )" +fc T ^{fi, oj)zW L , k ((3, w) 

n— 1 n— 1 fe=l 

+ R% N iP,w)zW L (P,u). (3.31) 

The last two lines will give a remainder i?^ 4< ^(/3, w), whose properties can be read out of those of 

the previous ones. This remainder has the same properties as R^ n (0,u>), and in fact it is a part 
of it. The rest of the proof is just algebra, and amounts to identify the right factors which enter 
the definition of QL,,rSfii w ) an d the expression of the full remainder. Here one must use (|3.14p and 
the notation introduced in (|3.15|) . (|3.23p and (|3.24|) . The proof is over. □ 
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4 Expansion of the trace of g^. The uniform bound. 



We have almost all ingredients needed for proving ()1.15p . Let (3 > 0, u>q > and to g as in the 
Section 1.1. We now need to take the trace in (|3.27p . Let us begin with the trace of the operator 
g L}0 ((3,u) (see (PHgj) ). If we use ([3~H?)) . and then ([3T4]) , ([3~TTj) and (pT2"3"|) , we can write: 

<?z,o(M = ff L (/3,w) - C^UM - r X K - ^(A^r 1 ^^ wJfiGS, w) 

iV 
fc=0 

~ C^R^UP, w)W L {fi, w)fxC9, w). (4.1) 



Apriori, this identity only holds in the bounded operators sense. But we know that Wl(/3, oj) is a 
trace class operator. It means that the operator 

M(/?,w) :=g L (/3^)-r 1 fL(P,^) (4.2) 



is a trace class operator, since all other operators in (|4.ip are trace class. Note that the two 
individual terms in M(/3, u>) might not be trace class. Now since M(j3, u>) has a continuous integral 
kernel M(-, (see Lemma |3~T1 and (|3.8[> ). its trace will be given by: 



Tr M(j3,w)= / M(x,x;/3,w)dx (4.3) 
Ja 

= / g L fax;P,u)dx- f -1 / f L (x, x; /3; w)dx 

Ja 

= / 5L(x,x;/3,w )<ix - / ri(x,x;/3,u;)dx. 

JA JA 

The last line is very important, since it shows that the "tilde" disappears when we take the 
trace. This is because the magnetic phase <j>(x, x) = for all x. But now J A <7l(x, x; 0, ojo)dx = 
Tr gL(f3,u>o), and we here recognize the very first term on the right hand side of (|1.13j) . Now if 
we use (I3.6[) p.8[) and (|3.9p we can write: 

N 

TtM(fi,u) =Tr 5i (/3,w )-r i y](M n / r L ,„(x,x; /3,w )dx 

n =i ^ 

+ (M W+1 < 1) (/5^,A r ), (4.4) 

where lu -> n^\/3,uj,N) is a smooth function in w near wo- Moreover, due to (|3.13p we obtain 
that the above integrals grow at most like -L 3 , as required. 

But there are several other terms which remain to be considered in (|4.ip and (|3.27|) . They are 

respectively Tr {S^kifi, u)Wl(0, u)r L (u)} and Tr g L ,n(0, w). 

These traces have two important things in common. First, we always take the trace of a 
product of integral operators with continuous kernels. Second, they all still depend on 5u>, but 
only through the magnetic phases; all factors are regularized operators, as defined in (|2.2p . We 
will now try to discuss all these different terms in a unified manner. 

Fix u>o > 0. Let wet and 5uj as above. Consider a product of operators of the form 

T(oj) := T (uj)Tt(w)...T n (w) 

where Ti{uj) are the regularized operators associated to some integral operators Tj(u/o), i = 0, n 
(see (|2.2|l ) and assume that this product is of trace class. Denote by t,(-, •) the kernel of Ti(wo), 
which is supposed to be jointly continuous in x and x'. We denote by fl„ the following flux related 
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quantity 

n-1 

Ai(x,yi) = 0, fl„(x,yi, ...,y„) = ^ fl(x, y k , y k+1 ) 

k=l 

= </ , (x,yi) + 0(yi,y 2 ) + --- + 0(yr l -i,yn) + 0(yT l ,x), n > 2. (4.5) 

Another important property of these operators is that their kernels are exponentially localized 
near the diagonal (see ()2.5f> . (|3.ip and (|3.13p ). Therefore there exist two positive constants C and 
a, independent of L, such that: 

max\U(x,x')\ < Ce- Q l x - X 'l, (x,x') e A 2 . (4.6) 
Then the diagonal value of the kernel of T{ui) reads as 

T(x,x,w)= / d yi ... / dy„e^ fl "^ 1 '-^^ (x,y 1 , Wo )i 1 (y 1 ,y 2 ,^ )... 

w )i„(y„,x, w ), (4.7) 

where we added together all individual phases from each regularized factor. Because we assumed 
that T is a trace class operator, the trace of T is 

TrT(w)= / T(x,x,cj)dx. (4.8) 

J A 

For m > 0, n > 1, let us introduce the notation: 

d m , n {L) := dx dyi... / dy„[ifi n (x, y lf y„)] m t (x, y 1 , uj ) 
Ja J a J a 

*i(yi,y2, w )...i„_i(y„-i, y„, w )t„(y„, x, uj q ). (4.9) 

Lemma 4.1. For every m > and n. > 1, t/iere exists a constant independent of L but depending 
on m, n such that 

\d m , n {L)\ < consti 3 . (4.10) 

Moreover, for a given N > 1 we have 

N 

Tr T(w) = ^ (^) m d m ,„(i) + (<5a;) JV + 1 7e i (a;, AT), (4.11) 

m=0 

where u> — > Hl(uj, N) is a smooth function near ujq. 
Proof. The equality (14. lip comes straight out of (|4.7p and 



Now let us prove the estimate (|4.10p . We recall the following estimate (|3.7p on the magnetic 
flux, 

|fl(x,y,z)| < |x-y||y-z|. (4.12) 
Then by induction one has for all n > 1, 

|fl„(x, yi , ...,y„)| < (|x - yi| + |yi - y 2 | + ... + |y n _ a - y„|) 2 . (4.13) 

Therefore the polynomial growth induced by this flux is diagonalized, i.e. it only depends on 
differences between the variables yi, j/i+i- But due to l|4.6p . we can write 

\d m , n (L)\ < const dx dy ± ... / dy„e~^ |x ~ yi1 
Ja ja ja 

e -f |yi-ya| e -f |yn-i-yn| e -f lyn-x^ (4-14) 

for some L independent constant. Here we used the exponential decay to bound the polynomial 
factors and then we extend the y integrals to the whole E 3 . So the volume growth is only given 
by the integral over x in the r.h.s of (|4.14p . The proof the lemma is over. □ 
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4.1 Proof of ( TTT3I) and of Theorem Q 

We can now put together the results of this section and prove the key estimate (|1.13[> . In Theorem 
13.51 we obtained an expansion for <7z,(/3, u>) as announced in (|1.14p . When we take the trace of 
gi,(f3,<jj), the term R^ N (/3,u>) will only give a contribution to the remainder in (|1.13p . hence we 
ignore it. 

Then the term gL,o((3, <-u) given in (|3.28|) can be written as a sum between an operator M(f3, u) 
from ()4.2|) . and a sum of operators of the type treated in Lemma \A. II Then from (|4.4p and the 
above mentioned lemma we can conclude that (d^Tr <7l,o)(/3, <^o) grows at most like L? . 

Finally, looking at the contribution coming from QL.n 

((3,u>), with n > 1. Using the same 
lemma, we obtain in a similar way that (d^Tr gL,n)(P,wo) grows at most like the volume. We 
therefore conclude that (c^Tr gi){[},ujQ) behaves like L 3 , uniformly in £ and z, and the proof of 
(fl~T2|) is done. 
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